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ABSTRACT 

We report the detection of a narrow Fe Ka emission line in the Seyfert 1 
galaxy NGC 5548 with the Chandra High-Energy Transmission Grating. In 
the galaxy frame we measure a center energy of 6. 402to!o25 keV, a FWHM 
of 4515±2645 km/s, an intensity of 3.6l};5 x 1(T 5 photons cm~ 2 s~\ with an 
equivalent width of 133^54 eV (errors are 90% confidence for one parameter). 
The line is only marginally resolved at the 90% confidence level. The line 
energy is consistent with an origin in cold, neutral matter, but ionization 
states up to ~ Fe XVIII are not ruled out. We cannot constrain the detailed 
dynamics but assuming Keplerian motion, the velocity width is consistent with 
the line being produced in the outer optical/UV broad-line region (BLR) at 
about a light-month from the central X-ray source. We cannot rule out some 
contribution to the narrow Fe Ka line from a putative, parsec-scale obscuring 
torus which is postulated to be a key component of AGN unification models. 
The continuum intensity during the Chandra observation was a factor ~ 2 less 
than typical historical levels. If the X-ray continuum was at least a factor of 
2 higher in the recent past before the Chandra observation and the narrow 
Fe Ka intensity had not yet responded to such a change then the predicted 
line intensity and equivalent width for an origin in the BLR is within the 
90% measurement errors. Anisotropic X-ray continuum illumination of the 
BLR and/or additional line emission from a torus structure would improve the 
agreement with observation. Two out of three archival ASCA data sets are 
consistent with the narrow line being present with the same intensity as in the 
Chandra observation. However, there is strong evidence that the narrow-line 
intensity varied and was unusually low during one of the ASCA campaigns. In 
general, inclusion of the narrow line to model the overall broad Fe Ka line profile 
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in terms of a rotating disk plus black-hole model can have a non-negligible effect 
on the disk-line intensity and variability properties. Variability of the broad disk 
line in NGC 5548 is difficult to reconcile with the expectations of the simple 
disk model, even when the narrow-line component is accounted for. It will be 
important to ascertain the importance of a similar non-disk Fe Ka line in other 
Seyfert 1 galaxies. Future monitoring of the narrow Fe Ka component with large 
collecting area and high spectral resolution will enable reverberation mapping 
of the BLR region complementary to similar studies using the optical/UV lines 
and therefore provide independent constraints on the black-hole mass. 

Subject headings: galaxies: active - galaxies: emission lines - galaxies: 
individual: NGC 5548 - X-rays: galaxies 
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1. 



INTRODUCTION 



The Advanced Satellite for Cosmology and Astrophysics (ASCA - Tanaka et al. 1994) 
found that the Fe Ka fluorescent emission line in Seyfert 1 galaxies is often very broad and 
is generally interpreted as the result of an origin in matter in an accretion disk rotating 
around a central black hole (see Fabian et al. 2000 and references therein). The line 
profile is sculpted by characteristic gravitational and Doppler energy shifts. Currently, 
study of the Fe Ka emission line is the only way to probe matter within a few to tens of 
gravitational radii of a black hole. There is an important caveat with respect to studying 
the relativistic Fe Ka line profiles with ASCA. The energy resolution of the CCD detectors 
is only E/AE ~ 40 (or FWHM ~ 160 eV, or ~ 7500 km/s) in the Fe Ka band (6-7 keV). 
Thus, even for the highest signal-to-noise data, one could not unambiguously separate that 
part of the line which originates in the accretion disk, from any other component originating 
in matter beyond the accretion disk, further away (e.g. see Weaver and Reynolds 1998; 
Done, Madejski, and Zycki 2000). Such line components will be much narrower than the 
disk line, due to the lower Keplerian velocities and small gravitational shifts. Indeed, such 
a narrow Fe Ka emission line from the putative obscuring torus is thought to be the key for 
unifying type 1 and type 2 AGN, and is expected from Seyfert 1 galaxies (e.g. see Krolik, 
Madau, and Zycki 1994; Ghisellini, Haardt, and Matt 1994). The torus is extended on 
the parsec-scale and would lead to a FWHM for the narrow Fe Ka line from this region 



of ~ 760y (M 8 /r pc ) km/s (assuming a virial relation and an r.m.s. velocity dispersion of 
V / 31^whm/2; see e.g. Netzer 1990). On the same scale, narrow-line emission may also 
originate in the extended warm scattering region which is also a feature of unified AGN 
models (e.g. see Krolik and Kallman 1987). A broader narrow-line contribution (1000s of 
km/s) might also be expected from the optical/UV broad-line region (BLR). 

As we shall demonstrate in this paper, accurate measurements of the narrow Fe Ka 
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line energy, width and intensity can provide a suite of powerful physical diagnostics of 
the central engine in AGN. Up to now, the expected narrow Fe Ka line has never been 
unambiguously detected in a Seyfert 1 due to limited instrumental energy resolution and 
sensitivity. However, strong evidence for its existence and physical parameters of the 
putative line have been deduced for several sources (e.g. Weaver et al. 1993; Yaqoob & 
Weaver 1996; Yaqoob et al. 1995,1996; Sulentic et al. 1998; Weaver & Reynolds 1998; 
Done et al. 2000). The narrow line appears to be most prominent in NELG such as NGC 
2992 (Weaver et al. 1996), which may sometimes have an underlying broad component 
(e.g. MCG —5-23-16, Weaver et al. 1997). The above evidence for the narrow Fe Ka line 
in Seyfert 1 galaxies has been tempered by the case of MCG —6—30—15 in which, during 
a flare, the entire broad Fe Ka line was highly redshifted, leaving no clear emission peak 
around 6.4 keV (Iwasawa et al. 1999). This has led to some debate about the existence of a 
narrow component. 

In this paper we report the detection of a narrow Fe Ka line in the Seyfert 1 galaxy 
NGC 5548 from a Chandra observation. In §||] we present the data and describe the analysis 
procedures. In §|] we show that the Fe Ka line in NGC 5548 cannot be explained by 
emission from a relativistic disk alone and present the basic measurements of the narrow 
Fe Ka line parameters. In §^ we discuss all the different possible mechanisms contributing 
to the width of the narrow line and give our best estimate of its velocity width. We discuss 
the physical constraints imposed by the narrow line in §|5], and in we investigate the 
effect of the narrow line on the parameters derived for the relativistic broad line using 
archival ASCA data. In §|7] we discuss some implications of the presence of a non-disk Fe Ka 
component on the interpretation of variability studies of the line and continuum. Finally, 
we present our conclusions in §|| 
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2. THE CHANDRA DATA 

NGC 5548 was observed with Chandra starting on 2000 February 5, UT 15:38 for 
a total duration of 83 ks (continuous). The instrument used in the focal plane of the 
High Resolution Mirror Assembly (HRMA) was the High-Energy Transmission Grating 
(or HETG - Markert, et al. 1995). The HETG consists of two grating assemblies, a 
High- Energy Grating (HEG) and a Medium- Energy Grating (MEG). The HEG has a 
better efficiency than the MEG above ~ 5 keV. The spectral resolution of the HEG is 
roughly two times better than the MEG (AA = 0.012A and 0.023A for the HEG and MEG 
respectively). AA is roughly constant over the instrument bandpass (HEG: ~ 0.8 — 9 keV , 
MEG: ~ 0.5 — 9 keV). Events dispersed by the gratings are collected by a CCD array and 
can be assigned an energy based on the position along the dispersion axis. Since the CCDs 
have intrinsic energy resolution, background events can be rejected with a high efficiency, 
and different spectral orders can be easily discriminated. 

Our Chandra observation of NGC 5548 yielded a net exposure time of 82,198.4 s. After 
accounting for a deadtime factor of 0.987337, the effective exposure time is 81157.6 s For 
our analysis we used data from events reprocessed using the latest software and calibration, 
as of 3 July 2000 []. Genuine photon events collected by the CCDs fall into specific pixel 
patterns classified by their grade and we retained only grades 0, 2, 3, 4, and 6; rejected 
grades correspond to non-X-ray events. We examined the radial profile of the zeroth-order 
image and found a single Gaussian to be a poor representation, due to the effects of pile-up 
broadening the profile (the best-fitting width was a = 0.52 arcseconds). A double Gaussian 
model gave a good fit and the best-fitting widths of two Gaussian were 0.48 ± 0.01 and 
1.09 ± 0.07 arcseconds, the ratio of the normalization of the narrow to broad Gaussian 

4 Details of the Chandra pipeline reprocessing and data analysis specifics can be found at 
http: / / asc.harvard.edu/ciao/thrcads/ 
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being 27.4. Since the zeroth order data suffer from pile-up we did not use these data for 
spectral analysis. For the grating data, in the cross-dispersion direction we found that 
the counts distribution as a function of cross-dispersion angle could be approximated by 
Gaussians with widths in the range 0.47-0.58 arcseconds for the four arms of first-order 
data (±1 orders for the HEG and MEG). We accumulated grating spectra from events along 
the dispersion direction, and within ~ ±3.6 arcseconds of the peak in the cross-dispersion 
direction. In the first-order spectra we obtained total count rates (summed over the ±1 
orders) of 0.19 cts/s and 0.43 cts/s for the HEG and MEG respectively. The count rates in 
the higher order spectra are at most 11.3% and 9.8% of the first order count rates for HEG 
and MEG respectively so in this paper we will restrict our analysis to the first order spectra 
only. 

We do not expect further refinements to the CCD gain to affect the energy scale 
of the grating data since the energy assignment of events depends only on the position 
of the event on the CCD according to the dispersion relation. The CCD energy scale 
may, however affect the rejection of source events as background events and vice versa. 
Since we are most concerned with the Fe Ka band and this band will have the highest 
proportion of background in the data that is relevant for this paper, we checked how 
important background is for our data in the 5.5-7 keV band. We collapsed the 5.5-7 keV 
first-order HEG data in the dispersion direction and examined the count distribution as a 
function of cross-dispersion angle, which is again approximately Gaussian. We estimated 
the background level by fitting a model consisting of a Gaussian plus a constant and found 
that the background is ~ 1.5% of the total counts in the 5.5-7 keV band. The same 
procedure for the MEG showed that the background is ~ 2.2% of the total counts in the 
5.5-7 keV band. 

We made effective area files (ARF, or ancillary response file) using ciao v 1.1 (Chandra 
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Interactive Analysis of Observations), which also takes account of the dithering of the 
satellite. We only used CIAO for the reprocessing of the event data and to make the 
ancillary response files; the remaining data analysis (except for spectral fitting) was done 
using IDL routines written by ourselves. For the spectral fitting procedures described 
below, we used spectral response matrices acisheglD1999-07-22rmf N0002 . f its (HEG) 
and acismeglD1999-07-22rmfN0002.f its (MEG). Note that absolute flux- calibration may 
be uncertain by as much as 30% at the time of writing (Marshall et al. 2000). However our 
principal results (i.e. center energy width and equivalent width of the Fe Ka line) do not 
depend at all on the absolute flux calibration. 

The spectral energy resolution is determined principally by the telescope response 
function but could be degraded by grating period variation (only at very low energies), 
incoherent dispersion (~ 1% of first order events), incorrect aspect solution, and defocussing. 
By examining the zeroth-order image from CCD frame-transfer events (which do not suffer 
pile-up), Paerels et al. (2000) determined from a HETG observation of Cyg X-3 that these 
effects do not compromise the spectral resolution within the measurement uncertainties. 
We also examined the frame-transfer events (or 'streak data') of our NGC 5548 data and 
found that the photon distribution in the dispersion direction could be approximated well 
with a single Gaussian, resulting in a best-fitting width of a = 0.41 ± 0.03 arcseconds (90% 
confidence errors), consistent with nominal expectations. This reassures us that there are 
no adverse effects on the spectral resolution of the gratings. The absolute calibration of 
the wavelength scale is typically ~ 2mA (e.g. Paerels et al. 2000), much smaller than the 
statistical uncertainty in the Fe Ka emission line energy that we are concerned with in this 
paper. 
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3. DETECTION OF A NARROW Fe Ka LINE IN NGC 5548 

We used XSPEC for spectral fitting the 2-7 keV HEG and MEG spectra (the 
signal-to- noise for both the HEG and MEG drops rapidly above ~ 7 keV). The —1 and +1 
orders were summed for the HEG and MEG. The spectral bins had a width of ~ 0.005A 
which corresponds to the la width of the spectral resolution of the HEG. For most of the 
2-7 keV bandpass, in the HEG or MEG spectra, there were less than 10 counts per bin and 
at most 14 or 12 counts per bin (HEG and MEG respectively). Since the data are in the 
Poisson regime we performed the spectral fitting and error analysis using the C-statistic. 
Figure la shows the HEG and MEG spectra fitted with a simple power-law model, which 
describes well the overall spectra in the 2-7 keV band. The consistency between the HEG 
and MEG is good over most of this energy range, except that between ~ 4 — 5 keV the 
MEG shows systematic deviations of the order of ~ 20 — 30%. Figure la also shows that 
an Fe Ka emission line is clearly detected with high signal-to-noise in both the HEG and 
the MEG. The insets in Figure 1 show close-ups of the Fe Ka region for the HEG data, 
and also the Fe Ka region from an ASCA observation in July 1996. The ASCA Fe Ka 
line is broad, with no obvious signature of a narrow core (but note that the ASCA energy 
resolution in the Fe-K region is ~ 160 eV FWHM, a factor of ~ 4 worse than the HEG data 
in that band). 

When we added a Gaussian to the power law continuum to model the Fe Ka emission 
line in the HEG plus MEG spectra, the C-statistic decreased by 25.6, which corresponds to 
a detection significance of nearly 5a for three degrees of freedom. In this latter model, there 
were six free parameters in total; two independent normalizations for the continuum (one 
for the HEG and one for the MEG), the power-law photon index (r), the line intensity, 
Gaussian width and center energy. The results of this spectral fit are shown in Table 1, 
along with the 90% confidence, one-parameter errors for the derived model parameters. 
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Also shown in Table 1 for reference are the spectral fitting results obtained when the same 
model is applied to the HEG data only. These latter results are consistent with the joint 
fit so hereafter we will refer only to the parameters derived from the combined HEG and 
MEG fits. We obtain a photon index for the power law of 1.54 ± 0.06 and a 2-10 keV flux 
of 2.4 x 10~ n ergs cm~ 2 s _1 . Historical 2-10 keV flux levels for NGC 5548 measured by 
EXOSAT, Ginga, and ASCA, according to the archival HEASARC databases and published 
results, span the range 1.7 — 6.4 x 10~ n ergs cm~ 2 s -1 , with values in the upper third of 
this range being most common (e.g. see Turner & Pounds 1989; Nandra & Pounds 1994; 
Chiang et al. 2000). Thus it appears that the 2-10 keV flux during the Chandra observation 
was a factor of ~ 2 down from typical historical values. 

Unless otherwise stated, we will refer all our measurements for X-ray line parameters to 
the galaxy frame, using a redshift of 0.01676, obtained from the NASA/IPAC Extragalactic 
Database (NED). Note that the redshift of NGC 5548 is uncertain by the equivalent of 
~ 100 km/s, and the value z = 0.01676 refers to the value otained from optical emission-line 
measurements, whilst Hi measurements give z = 0.017175. See, for example, Mathur, Elvis, 
and Wilkes (1999) and references therein. Thus, our quoted galaxy-frame Fe Ka line-energy 
measurements carry additional uncertainties of the order of ~ 2 eV. 

Table 1 shows that for the Gaussian fit we obtain a peak energy, intensity, and 
width of the narrow Fe Ka line detected by the HEG and MEG, of 6.402^;^ keV, 
3.6+^5 x 10~ 5 photons cm~ 2 s _1 , and 4515+2645 km/s FWHM respectively. The equivalent 
width (EW) of the line is 133tg eV. Statistical errors are 90% confidence for one parameter 
for all the above measurements. Yaqoob (1998) noted that, in the extreme Poisson limit, 
the C-statistic can under-estimate the statistical errors on spectral features which are 
comparable to or narrower than the instrument resolution. Although our spectra are not 
in the extreme Poisson limit, we performed simulations (as described in Yaqoob, 1998) to 
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check the reliability of the statistical errors derived from spectral fitting. Thus, we made a 
grid of simulated spectra in which the input Fe Ka line intensity was set at nine values in 
the range 1.5-6.0 xlCT 5 photons cm~ 2 s^ 1 and the input (Gaussian) line width was set at 
twelve values in the range 4-80 eV. The input line energy was set at 6.402 keV in the galaxy 
frame (the best-fit value from spectral fitting). Each grid point contained 1000 simulated 
spectra. All other model parameters (e.g. exposure time) were set at the values obtained 
from the Chandra data. The simulated spectra were fitted with the same power law plus 
Gaussian model and we found the statistical errors on the derived model parameters to be 
consistent with those obtained from spectral fitting. Our simulations also show that the 
HEG data contain a mean of only 70 photons in the Fe Ka line with a la dispersion of ±8 
photons. This is consistent with the mean number of line photons (72) we obtain from the 
product of the effective area at the observed line center energy, the exposure time and the 
best-fitting line intensity from Table 1. A similar calculation shows that the expectation 
value of the number of line photons detected by the MEG is 23. 

We next addressed the question of whether the narrow Fe Ka line detected by Chandra 
is resolved. Note that the spectral resolution of the HEG at the observed Fe Ka line energy 
corresponds to ~ 38 eV FWHM, or ~ 1800 km/s FWHM. Starting with the best-fitting 
power law plus Gaussian model to the HEG and MEG data we constructed confidence 
contours using the C-statistic for the Fe Ka line intensity versus the line width. Figure 2 
shows the 68%, 90% and 99% joint confidence regions, from which it can be seen that the 
emission line is resolved at the 90% level but not at 99%. 

3.1. RELATIVISTIC DISK LINE CONTRIBUTION 

We note that the Fe Ka emission line detected by Chandra cannot be accounted for 
by a relativistic disk line alone. We fitted the HEG plus MEG data with a power-law 
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continuum and a model of the broad Fe Ka line originating in an accretion disk rotating 
around a Schwarzschild black hole (e.g. Fabian et al. 1989). We fixed the disk power-law 
emissivity index at —2.5 (close to the mean of a sample of Seyfert 1 galaxies - Nandra et al. 
1997), and the energy of the line in the disk rest-frame at 6.4 keV. The inner and outer disk 
radii were fixed at 6 rg and 1000r 9 respectively, where r g = GM/c 2 , and 6 rg is the radius 
of the last marginally stable orbit around a black hole of mass M. The data are not very 
sensitive to the parameters that were fixed. The inclination angle (0 o bs) and the disk-line 
intensity were allowed to float. Including two independent continuum normalizations for 
the HEG and MEG, there are only a total of five free parameters, one less than the power 
law plus Gaussian line model. The best fit was much worse than the one in which the 
Fe Ka line was modeled only by a Gaussian (see Table 1); specifically the C-statistic was 
higher by 10.1 for the disk-line model, for one less degree of freedom. This means that the 
confidence level for the Gaussian model being a better description of the line is greater 
than 99.8%. The best-fitting disk inclination was 14.6^4;q degrees (90% confidence for one 
parameter), which is inconsistent with the best-constrained inclination derived from fitting 
ASCA data (see §|6|). Essentially, the narrow Fe Ka line detected by Chandra is too sharply 
peaked for even low-inclination disk-line models. 

Next, we investigated the effect of including a broad, relativistic Fe Ka line in addition 
to the narrow Gaussian model. This time, in addition to the parameters which were fixed in 
the disk-line only fits, we also fixed the disk inclination angle at the best-constrained value 
from the ASCA data, of 31° (see All the Gaussian parameters were still free so there is 
only one additional degree of freedom in this composite, dual-line model. Spectral fitting 
the HEG and MEG data with this model reduces the C-statistic by only 0.1 relative to the 
power law plus Gaussian model. We note here that the Chandra data are insensitive to the 
disk inclination angle, with the C-statistic changing by less than 0.6 when the inclination 
angle is varied over the entire range of 0° to 90°. Since the relativistic line is so broad, its 
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flux is spread out over a wide energy range and we found that it has a negligible effect 
on the derived parameters of the narrow Fe Ka line. In particular, the intensity of the 
narrow line decreased by ~ 3% and the width increased by ~ 4% when the disk line was 
included. The best-fitting intensity and equivalent width of the relativistic disk line was 
0.4i£i x 1(T 5 photons cm 2 s 1 and 13-il 7 e V respectively (errors are 90% confidence for 
three parameters). Thus, the broad component of the Fe Ka line is not detected due to the 
small effective area of the Chandra gratings but the statistical upper limits on its intensity 
are consistent with historical ASCA measurements (Chiang et al. 2000, and §|6]). 

4. BROADENING MECHANISMS FOR THE NARROW Fe Ka LINE 

The intrinsic width of the Fe Ka line, which results from a combination of the natural 
width and level width, is less than 1 eV (Bambynek 1972) so for our purpose it is negligible. 
The thermal width of the narrow Fe Ka line is ~ 1.7-^/ (T/IO^K) km/s so this too is 
negligible. We can also neglect gravitational redshifts (these are of order AE/E ~ r g /r for 
r/r g <C 1, where r g = GM/c 2 ) which will be less than 1 eV at a distance greater than four 
light days from a black hole with a mass of 1O 8 M ). In §|5| and Figure 3 we show that the 
narrow Fe Ka line originates at ~ 1 light month from the central X-ray source. 

The Fe Ka line actually consists of two lines, Ka\ at 6.404 keV and Ka 2 at 6.391 keV, 
with a branching ratio of 2:1. Therefore we investigated how this affects our results as 
follows. Using the 2-7 keV HEG and MEG data we fitted a model consisting of a power-law 
continuum plus two Gaussians, with the center energies fixed at the above energies, and 
the intensity of the Gaussian at higher energy was constrained to have twice the intensity 
of the other Gaussian. The widths of the two Gaussians were constrained to vary together. 
We obtained best-fitting values for the total line intensity (the sum of the two Gaussian 
intensities) and the width which were the same as those obtained from a single Gaussian 
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model (Table 1), and also obtained similar statistical errors. In particular we obtained a 
width of the line of 4515 j^g km/s FWHM (compare with the single- Gaussian measurement 
in Table 1). 

If the narrow Fe Ka line is formed in Compton-thick material such as the putative 
obscuring torus, Compton down-scattering will also increase the apparent width of the line. 
However the Compton shoulder on the red side of the line peak due to the first scattering 
is already an order of magnitude down in intensity relative to the unscattered component 
(e.g. Matt, Brandt, and Fabian 1996; Iwasawa, Fabian, and Matt, 1997). The present data 
are not sensitive to this level of structure. 

Another uncertainty arises from a possible contribution to the Fe Ka line from higher 
ionization states such as Fe II, which is expected if some of the line is coming from the 
putative obscuring torus. However, laboratory measurements of the energies of these higher 
ionization lines is uncertain and the HEG does not have sufficient energy resolution to place 
any tighter constraints other than the fact that the 90% confidence errors on the line center 
energy rule out a dominant contribution from ionization states higher than Fe XVIII (e.g. 
Makishima 1986). Thus, until better data are available we attribute the broadening of the 
narrow Fe Ka line in the Chandra data for NGC 5548 to the bulk velocity of the emitting 
gas but do not rule out some contribution from emission components with velocity widths 
less than our formal lower limits imply. 



5. PHYSICAL DIAGNOSTICS FROM THE NARROW Fe Ka LINE 



The narrow Fe Ka line discovered in NGC 5548 immediately enables some robust 
deductions to be made. We showed in § [3.1| that line is too narrow to be the core of a 
relativistic disk line (e.g. see Fabian et al. 1989, 2000) in which the line-emissivity falls 
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off with increasing distance from the continuum source. In the following sections we use 
the measurements of the peak energy, width, intensity and equivalent width of the narrow 
Fe Ka line to place physical constraints on the line-emission region. 

5.1. PEAK ENERGY OF THE LINE 

The peak energy of the line is ~ 6.40 ± 0.03 keV, consistent with the rest-energy for 
cold, neutral iron. The upper limit on the energy corresponds to an ionization state of 
Fe XVIII and this is consistent with the ionization parameter for the optical/UV BLR 
deduced by Goad and Koratkar (1998) of \ogU ~ —0.6 — 0.0 from detailed modeling. As 
mentioned in §^j, gravitational redshifts are not large enough to significantly shift back 
the line energy from ionized Fe to lower energies (the distance of the line-emitting region 
we deduce below produces energy shifts less than 1 eV). Therefore a dominant emission 
component from a warm scattering zone, consisting of He-like or H-like Fe, (a component 
of type 1 and type 2 AGN unification models, e.g. Krolik & Kallman, 1987) appears to be 
ruled out, at least in NGC 5548. Models in which the broad Fe Ka line is produced and 
broadened by Compton scattering in an optically thick medium (e.g. Misra and Kembhavi 
1998; Mouchet et al. 2000) are likewise ruled out since the unscattered peak of the Fe Ka 
line would be at ~ 6.67 keV. 

We can also rule out a scenario in which the narrow Fe Ka line is not seen directly but 
from reflection in optically thin cold gas out of the line-of-sight. Line photons of energy Eq 
will be Compton down-scattered to a mean energy of Eq/[1 + Eo/m e c 2 }. For Eq = 6.4 keV 
this is ~ 80 eV lower than the peak energy we measure from the HEG data for NGC 5548. 
Of course if the scattering medium is hot, Fe XX could produce a line which is scattered 
back down to the observed energy. However, this requires a temperature of the order of 1 
keV, and this would give rise to a Compton-broadened line with a FWHM ~ 670 eV, an 
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order of magnitude larger than observed (see Pozdnyakopv, Sobol and Sunyaev, 1983). In 
fact, to obtain a line width consistent with the observed value, the gas electron temperature 
would have to be less than 1.4 x 10 5 K. 

5.2. LINE WIDTH 

From a detailed study of the velocity widths and line/ continuum time lags for 
optical/UV emission lines in NGC 5548, Peterson & Wandel (1999) found strong evidence 
for Keplerian motion of the line-emitting matter and estimated the central mass to be 
6.8 x 10 7 M Q . As Peterson & Wandel (1999) caution, this mass estimate is subject to 
uncertainties in the geometry and kinematics of the BLR, and optical/UV line reprocessing 
physics. The actual mass may be as small as ~ 5 x 10 6 M and as high as ~ 10 8 M Q . A 
summary of estimates of the central mass in NGC 5548 using other methods can be found in 
Hayashida et al. (1998) and are consistent with this general range. In Figure 3 we plot the 
Fe Ka line width constraints (90% confidence range) that we measured from the Chandra 
data on a diagram of reverberation time-delays versus velocity widths of some of the BLR 
optical/UV emission lines measured by Peterson & Wandel (1999). Using their value of 
the central mass and our Fe Ka FWHM of 451512645 km/s we obtain a size of the Fe Ka 
line-emitting region of 23 light days (90% confidence, one parameter limits: ~ 7 — 135 light 
days). This places the Fe Ka line-emitting region predominantly in the outer BLR and 
beyond (from Figure 3 the BLR extends from ~ 1 — 30 light days). Since we infer a low 
ionization state for the Fe Ka line-emitting region we would also expect Mg II A2798 from 
the same region. Indeed, Goad & Koratkar (1998) measure a FWHM of 358liif[ km/s 
for the broad core of Mg II A2798. It is not clear why the innermost region of the BLR is 
apparently contributing little to the Fe Ka line. Perhaps the BLR covering factor increases 
radially and is significantly less in the inner BLR than it is in the outer BLR. The deduced 
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location of the narrow Fe Ka line emission also implies that we might not expect the line 
intensity to respond to changes in the X-ray continuum intensity on timescales less than 
a week or so. Also, the energy resolution of the HEG does not allow us to rule out some 
additional contribution to the Fe Ka line from material further out, such as the putative 
obscuring torus. 

Yaqoob et al. (1993) calculated the expected line profiles, equivalent widths, and 
transfer functions for an Fe Ka line formed in the BLR. For pure radial motion in a 
spherically symmetric cloud distribution (inflow or outflow), the line profile has almost 
vertical sides, nearly the same intensity on the red and blue sides for the narrow widths we 
are considering here. The full width of such a line profile corresponds to twice the maximum 
flow velocity. We fitted the HEG plus MEG data with a very simplistic rectangular model 
of the line profile and derive a flow velocity of 2910ti4go km/s. 



5.3. LINE INTENSITY AND EQUIVALENT WIDTH 

We estimate the intensity and equivalent width of the Fe Ka line from the BLR 
following Krolik and Kallman (1987) but use updated Fe abundance and K-shell cross- 
section data. We assume a spherically symmetric cloud distribution and neutral iron 
for simplicity. The line intensity can be obtained by computing the number of photons 
removed from the continuum above the Fe-K edge threshold energy (at E K , or 7.11 keV 
for neutral Fe) and multiplying by the fluorescence yield (ujk = 0.34 for neutral Fe) and 
the fraction of the sky covered by BLR clouds (/ c ). If the observed power-law continuum is 
N P E~ T photons cm" 2 s _1 and ope-K^) is the K-shell photoelectric absorption cross-section 
as a function of energy, then the predicted observed intensity of the Fe Ka line is 



- 18 - 



roo 

I Fe - K = fc^Kha / N p E~ r [l - exp (-a Fe _ K A Fe N H )} photons cm -2 s" 1 (1) 

J E K 

where N H is the column density of each cloud, A Fe is the Fe abundance relative to 
Hydrogen, and /k is the fraction of emission- line photons appearing in the Fe Ka line (as 
opposed to Fe K(3 line), and we take /k = 150/167 (e.g. Bambynek 1972). Assuming the 
clouds are optically thin, we use a linear approximation for the exponential in equation 
1. We fitted photoelectric absorption cross-section data from Henke tables^ with a simple 
power-law and obtained 

/ R \ -2.646 

a Fc ^ K (E) = 3.79 I — I x 1(T 20 cm 2 per atom. (2) 
\E K J 

for E > Ek- This is not as steep as the (E/Ek)~ 3 relation used in older work (e.g. 
Krolik and Kallman 1987; Yaqoob et al. 1993). Using an Fe abundance of 4.68 x 10~ 5 
(Anders and Grevesse 1989) the line intensity is then 



T Fe _ Ka = 1.04 x 10 1 1 1 ] 1 ! 1 



0.35/ V0.34/ V4.68xl0- 5 / V10 23 cm 2 
12 X / Ek ^ T (7.11) L5 - r ( ACO N \ n _, ) photons cm" 2 s" 1 . (3) 



T + 1.646/ V7.ll/ V4.68 x 10 

The equivalent width of the line is 



EWi 



Fc-Kq 




( ^ 



H 



Vio 23 



cm' 



x 



(4) 



5 http: / /www-cxro. lbl.gov/optical_constants 
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Here E Ka is the center energy of the Fe Ka line. All quantities in the above equations 
refer to the source frame. The above approximations are valid if the Fe-K absorption optical 
depth at all energies is much less than unity. Effectively this means Nh <C 5.6 x 10 23 cm -2 
(using the maximum optical depth, at Ek)- Exact calculation using numerical integration 
showed that for N H = 10 23 cm~ 2 , the approximations lead to underestimates of the line 
intensity and equivalent width of ~ 6%. 

In a recent detailed study involving modeling the broad optical/UV emission lines 
in NGC 5548, Goad and Koratkar (1998) found that a covering factor of 35% and a 
column density of 10 23 cm -2 satisfy the observational constraints. These values, along with 
parameters measured for NGC 5548 from the Chandra data, namely N p = 4.68 x 1CT 3 , 
T = 1.54, give a predicted Fe Ka line intensity of ~ 1CT 5 photons cm -2 s _1 and an EW 
of 36 eV. These both fall short of the observed values (3.6±i; 5 photons cm~ 2 s _1 for the 
intensity and 133^|| for the EW) by a significant factor. If the BLR is responsible for the 
production of the narrow Fe Ka line, the observed intensity and equivalent width could be 
reconciled with the predicted values if the continuum level was a factor of ~ 3 higher in 
the recent past before the Chandra observation such that the line had not responded to the 
decline in the continuum by the time of the Chandra observation. Historical values of the 
2-10 keV flux range from ~ 1.7 - 6.4 x 10~ n ergs cm" 2 s" 1 (from HEASARC archives and 
Chiang et al. 2000 and references therein), so such a scenario is possible. There is of course 
the possibility that there is an additional contribution to the narrow Fe Ka line from the 
putative obscuring torus of a few tens of eV. A torus component of the line is not ruled out 
by our velocity-width measurements and would bring the observed and predicted strength 
of the line into better agreement with a more modest implied change in the continuum 
level. Another effect which would increase the predicted equivalent width of the Fe Ka line 
is anisotropic X-ray continuum illumination of the BLR (see Yaqoob et al. 1993 for details). 



-20- 



6. EFFECT OF THE NARROW LINE ON THE PARAMETERS OF THE 

RELATIVISTIC LINE 

NGC 5548, like many other Seyfert 1 galaxies, is known to have a relativistically 
broadened Fe Ka line (Figure 1 (b); Mushotzky et al. 1995; Nandra et al. 1997). The broad 
line cannot be clearly seen in the Chandra grating data because of the small collecting area 
and in §[| we placed formal limits on the intensity of the broad component from the Chandra 
data. We now investigate the effect on the relativistic line parameters derived from the 
ASCA data when a narrow Fe Ka emission-line component is included. There have been 
three observing campaigns of NGC 5548 with ASCA in 1993, 1996 and 1998. The latter two 
consisted of multiple observations spaced over one to three weeks and we used the summed 
data in each campaign. We fitted the 3-10 keV spectra with a power law plus a model for 
the Fe Ka emission line originating in a disk rotating about a Schwarzschild black hole as 
described in § [3.1| . No additional narrow Fe Ka line was included at first. As in § |3.1| we fixed 
the disk-line parameters to which the data are insensitive, namely the power-law emissivity 
index at —2.5, the inner and outer disk radii at 6 rg and 1000r 9 , respectively, and the energy 
of the line in the disk rest-frame at 6.4 keV. The inclination angle (0 o bs) and the disk-line 
intensity (Id) were allowed to float. Acceptable fits were obtained for all three spectra 
and the results are shown in Table 2. It can be seen that 8 b s is consistent with a single 
value (the weighted mean is 31.5° ± 6.5°) but the EW during the July 1998 observation 
was almost a factor 2 less than that in the other two epochs, implying time delays between 
variation in the continuum and line intensity. Indeed the line intensity was less in this 
observation than in the previous two observations which both had a lower continuum flux. 

Next, we added a narrow Gaussian to the model, with center energy and width fixed 
at the values derived from the HEG plus MEG measurements (Table 1). The intensities 
of both the broad and narrow components were allowed to float. The results are shown in 
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Table 3. The disk inclination angle is consistent with the spectral fits which did not include 
a narrow line (Table 2) but it is poorly constrained in the 1993 and 1998 observations. The 
best constraint, from the 1996 observation is o b s = 31 ± 8 degrees. It can be seen from 
Table 3 that all the ASCA data for 1993 and 1998 are consistent with the existence of a 
narrow Fe Ka line with the same intensity measured by Chandra. However, for the 1996 
observation, the 90% confidence, three-parameter errors on the narrow-line intensity do not 
formally overlap with the Chandra measurement. We conclude (with caution because of 
uncertainties in absolute calibration), that the narrow-line intensity is variable and was less 
during the 1996 ASCA observation than the values measured by Chandra and by ASCA in 
1993 and 1998. This conclusion is supported by the fact that A% 2 decreases by 8.1 and 76.9 
(for an additional degree of freedom) when a narrow line is included in the 1993 and 1998 
ASCA dual-line spectral fits but Ax 2 = for the 1996 ASCA data (see Table 3). 

7. IMPLICATIONS FOR VARIABILITY OF THE Fe Ka LINE IN 

SEYFERT 1 GALAXIES 

The intensity of the disk line derived from the ASCA dual-line spectral fits (§^) appears 
to be variable, it naturally being less for larger values of the intensity of the best-fitting 
narrow-line component (Table 3). Although it appears that both the intensity and EW of 
the disk line may be variable amongst the three ASCA observations, the 90% confidence 
errors formally overlap so we cannot draw any firm conclusions. If the disk-line variability 
is real then the line intensity is not simply correlated with the observed continuum level 
because the 1998 observation, the one with the highest continuum flux, appears to have 
the smallest disk-line intensity. Time delays between the continuum and line variability 
are expected to be less than 10 hours for a black hole mass of 7 x 1O 8 M (appropriate for 
NGC 5548; e.g. Peterson and Wandel, 1999) for line emission from within 100r 5 . Since 
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the continuum does not vary significantly on this timescale during the observing periods 
this implies that the response time of the disk line may be longer than expected, implying 
more emission from larger radii than is assumed. Thus, the broad disk-line variability, if 
it is real, has difficulty being reconciled with the simple disk model. Chiang et al. (2000) 
actually found that the total Fe-K line intensity (as measured by RXTE) is anticorrelated 
with the X-ray continuum flux. Now we have shown that when one takes account of the 
narrow Fe Ka line the problem is worse, since Table 2 shows that even without a narrow 
Fe Ka line the disk line EW in July 1998 is already less than expected ( and any additional 
narrow Fe Ka line can only decrease the disk-line intensity). 

Since the continuum level during the Chandra observation was at least a factor ~ 2 
less than that during the ASCA observations we might expect the disk line to have a 
smaller intensity in the Chandra data. However, in §|| we showed that the Chandra data 
cannot constrain the broad-line component very well and obtained an upper limit of 
7.6 x 10~ 5 photons cm" 2 s _1 on its intensity, formally consistent with the values obtained 
from all the ASCA observations. 

Recent campaigns to investigate the short-term variability of the Fe Ka line in Seyfert 1 
galaxies and its relation to the continuum variability have been conducted for a handful 
of sources (NGC 7314, Yaqoob et al. 1996; MCG -6-30-15, Iwasawa et al. 1996, 1999, 
Reynolds 2000; NGC 3516, Nandra et al. 1999; NGC 4051, Wang et al. 1999; NGC 5548, 
Chiang et al. 2000; NGC 7469, Nandra et al. 2000). If one assumes that all or most of 
the Fe Ka line emission comes from the putative accretion disk then one expects the line 
intensity to track the continuum on timescales appropriate for the size and geometry of of 
the X-ray source and disk system. Evidence for this is found only in NGC 7314 and NGC 
7469. The other sources show complex behavior which cannot be interpreted in terms of 
this simple model. Even in NGC 7314, only the broad wings of the Fe Ka line respond to 
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the continuum, leaving a constant intensity core centered at 6.4 keV (Yaqoob et al. 1996). 
This implies that NGC 7314 also has a separate, non-disk narrow Fe Ka emission line 
component, just like NGC 5548. The NGC 7469 data were taken with RXTE so do not have 
the spectral resolution to identify a constant core component of the emission line. Moreover, 
no study has yet shown the expected correlation between the Fe Ka line intensity and the 
Compton- reflection continuum (which also comes from the accretion disk). This is such a 
fundamental prediction of the disk model for the Fe Ka line, yet when it has been tested 
either no correlation is found, or else an anti- correlation has been found (such behavior is 
actually found in NGC 5548; see Chiang et al. 2000 and §^|). If a narrow, non-disk Fe Ka 
line component becomes important in general for low-flux states of Seyfert 1 galaxies, then 
variability studies need to take this into account. While discussions of a non-disk component 
to the Fe Ka lines have been around ever since their discovery in Seyfert 1 galaxies with 
early proportional counters in the late seventies, the Chandra observation of NGC 5548 now 
provides direct observational proof that this component exists, at least in one source. 



8. CONCLUSIONS 

Using the Chandra High-Energy Grating we detected a narrow Fe Ka emission line 
in the Seyfert 1 galaxy NGC 5548, with a center energy of 6.40 ± 0.03 keV, a FWHM of 
45151 2645 km/s and an equivalent width (EW) of 133^54 eV (90% confidence errors, three 
parameters of interest; all quantities are in the galaxy frame). The line is only marginally 
resolved. This emission line is distinct from the broad Fe Ka line in NGC 5548 and other 
Seyfert 1 galaxies, thought to originate in a disk rotating around a central black hole. 
The Chandra data are inconsistent with the line originating only from a disk because it is 
too sharply peaked even for low-inclination disks. The narrow Fe Ka line center energy, 
velocity width and equivalent width measurements are consistent with the bulk of the 
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line originating in cold, neutral matter in the outer BLR at about a light-week to a few 
light-months from the central source. We also do not, however, rule out ionization states 
up to about Fe XVIII. 

The agreement of the line equivalent width and intensity with predicted values 
assumes a spherically symmetric distribution of BLR clouds with a covering factor of 35% 
and typical cloud column densities of ~ 10 23 cm" 2 (values consistent with optical/UV 
emission-line studies) and is conditional upon the X-ray continuum having been a factor 
~ 3 higher in the immediate past prior to the Chandra observation. The 2-10 keV flux 
during the Chandra observation was a factor of 2.5 below typical historical values so this 
is not unreasonable. The continuum variability constraint would be relieved somewhat if 
there were a contribution to the narrow line from a putative, parsec-scale obscuring torus, 
which has been postulated as a key component of unification schemes of Seyfert 1 and 
Seyfert 2 galaxies. Indeed, such a contribution is not ruled out by the data. Anisotropic 
X-ray illumination of a spherical BLR would also increase the expected equivalent width of 
the narrow line. Considering its deduced location, the narrow Fe Ka line iteself may vary 
on timescales of weeks to months. Indeed, examination of archival ASCA data provides 
evidence for variability of the narrow Fe Ka line. 

We do not yet know whether the narrow Fe Ka emission line is a common feature of 
Seyfert 1 galaxies. From a moderate sample of sources, Nandra et al. (1997) placed an 
upper limit on a narrow Fe Ka component in the overall Fe Ka line profile typically of 
~ 60 eV, but some sources show stronger evidence for a separate, narrow Fe Ka line (e.g. 
Yaqoob & Weaver 1996; Weaver and Reynolds 1998; Done et al. 2000). In the case of NGC 
5548 one of the likely reasons for the relatively large EW of the narrow component observed 
by Chandra is that the continuum intensity happened to be significantly lower during the 
Chandra observation compared to the upper range of historical values. It will be important 
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to obtain Chandra data for a large sample of Seyfert 1 galaxies. However, studying the 
variability of the narrow line is rather expensive with Chandra . Monitoring the narrow 
Fe Ka line and using it for reverberation mapping with future X-ray astronomy missions 
with larger effective area will allow us to independently measure the central black-hole 
mass, complementing reverberation studies using optical/UV emission lines. Combining 
these future variability studies with spectral analysis with an energy of resolution of the 
order of 1 eV will allow us to constrain the details of the geometry and dynamics of the gas 
which are not accessible by current observations. 

We thank Dr. Kimberly Weaver and Professor Julian Krolik for valuable discussions 
relating to this work. We also thank Drs. H. Marshall and F. Nicastro for their helpful 
advice on Chandra data analysis issues. We are grateful to an anonymous referee for 
reviewing the manuscript and making some important suggestions. This research made use 
of HEASARC online databases. 



Table 1. Parameters of the narrow Fe Ka Line Derived from the HEG Data 



Parameter HEG only HEG and MEG 

Center Energy (keV) 6.390tffi 6.4021K 

Gaussian width (eV) 41±g 41^ 

Intensity (10~ 5 photons cm" 2 s" 1 ) 3.2±i;| 3.6±i;| 

Equivalent Width (EW) 1201$ 133lg 

Velocity FWHM (km/s) 4525±^gg 4525l|i 

Photon Index (r) 1.531±8;g|i 1.5401K 

C-statistic 842.7 1332.8 

degrees of freedom 1765 2648 

2-10 keV flux (10- 11 ergs cm" 2 s" 1 ) 2.5 2.4 



Note. — All line parameter values are referred to the rest frame of 
NGC 5548 (z = 0.01676). Statistical errors are 90% confidence for one 
interesting parameter (AC = 2.706). The 2-10 keV flux for the HEG 
plus MEG fit is the average from the two gratings. 
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Table 2. Relativistic Disk Line Fits to ASCA Data for NGC 5548 



Parameter July 1993 July 1996 July 1998 



Exposure (ks) 


29.6 


56.9 


130.0 


Disk inclination/ 9 ohs (degrees) 




3111 


29l? 2 


Disk-line intensity (10~ 5 photons cm" 2 s _1 ) 


10.3IL 7 


11.310 


7.6l|? 


Disk-line EW (eV) 


2341^ 


208±1° 


1201^ 


2-10 keV flux (KT 11 ergs cm" 2 s" 1 ) 


5.1 


5.4 


6.5 


x 2 


1016.3 


1199.8 


1483.9 


degrees of freedom 


1004 


1195 


1466 



Note. - 3-10 keV spectra fits using ASCA data (both SIS and GIS) for three 
observations of NGC 5548. The model consists of a power law plus an emission line (rest 
energy 6.4 keV) from a relativistic disk rotating around a Schwarzschild black hole. The 
inner and outer disk radii are fixed at 6 and 1000 gravitational radii respectively. The 
line emissivity is a power law with index fixed at —2.5. All the line parameters shown 
refer to the galaxy frame. All statistical errors on parameters are 90% confidence for three 
interesting parameters (A% 2 = 6.251). 
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Table 3. Dual Emission-Line Model: Relativistic Disk Line Plus Gaussian Fits to ASCA 

Data for NGC 5548 



Parameter July 1993 July 1996 July 1998 



Disk inclination/ 6> bs (degrees) 


4D -28 


3111 


C9+38 


Disk-line Intensity (1CT 5 photons cm~ 2 s _1 ) 


7 i +5.5 
'• ± -5.6 


n-3±2:g 


r. c+4.1 
0.0_ 4 2 


Disk-line EW (eV) 


147+114 

J-^' -116 


208±H 


82l| 


Narrow-Line Intensity (1CT 5 photons cm -2 s _1 ) 


q o+2.2 
°- z -2.8 


o.o±S:g 


2.9li;| 


Narrow-Line (eV) 


63li 


olo 7 


41+20 


x 2 


1008.2 


1199.8 


1407.7 


degrees of freedom 


1003 


1194 


1428 



Note. — Spectral fits to -AS'CA data for NGC 5548 using a composite model for the Fe Ka 
line, consisting of a relativistic disk line plus a narrow Gaussian. The Gaussian width of the 
narrow line is fixed at the best-fitting value obtained from the HEG plus MEG Chandra spectral 
fits (a = 41 eV, Table 1). All line parameters shown refer to the galaxy frame. All statistical 
errors are 90% confidence for three interesting parameters (A% 2 = 6.251). 
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Figure Captions 



Figure 1 

(a) Chandra HETG spectra for NGC 5548 from the HEG (solid error bars) and MEG 
(crosses for data points with dotted error bars). Solid line is the best-fitting power- law 
continuum to HEG plus MEG data but with the best-fitting HEG normalization. The 
±1 orders are summed. The narrow Fe Ka emission-line component is clearly detected in 
both the HEG and MEG. (b) ASCA SIS0+SIS1 relativistic Fe Ka line profile for the July 
1996 observation of NGC 5548 (see text), (c) Close-up of the HEG spectrum; the data 
are smoothed with a boxcar five bins wide, where the binsize is 0.005A. The solid line is 
the best-fitting model which consists of a power-law continuum and a narrow Gaussian. 
All parameters have the best-fitting values derived from fitting the HEG data (Table 1). 
Energy scale in all three panels is for the observer's frame. 

Figure 2 

Joint two-parameter confidence region of the narrow Fe Ka line intensity versus Gaussian 
width from spectral fitting to the HEG and MEG Chandra data. The contours correspond 
to 68%, 90% and 99% confidence. 

Figure 3 

Measurements by Peterson and Wandel (1999) of the time lags (relative to the appropriate 
continuum ) and the FWHM of various optical/UV emission lines in NGC 5548. The 
dashed line corresponds to the best Keplerian model fit obtained by Wandel and Peterson 
(1999) for the data from which they obtain a central mass of 6.8 x 10 7 M . The dotted 
lines show the best-fitting parameters and 90% confidence (one-parameter) limits on the 
FWHM of the narrow Fe Ka emission line measured from our Chandra HEG plus MEG 
data (see Table 1). 
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Fig. 1. — See figure caption. 
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Fig. 2. — See figure caption. 
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Fig. 3. — See figure caption. 



